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An angular effective mass formalism previously introduced is used to study the density of states 
in warped and non-warped energy bands. Band warping may or may not increase the density- 
of-states effective mass. Band “corrugation,” referring to energy dispersions that deviate “more 
severely” from being twice-differentiable at isolated critical points, may also vary independently of 
density-of-states effective masses and band warping parameters. We demonstrate these effects and 
the superiority of an angular effective mass treatment for valence band energy dispersions in cubic 
materials. We also provide some two-dimensional physical and mathematical examples that may 
be relevant to studies of band warping in heterostructures and surfaces. These examples may also 
be useful in clarifying the interplay between possible band warping and band non-parabolicity for 
non-degenerate conduction band minima in thermoelectric materials of corresponding interest. 
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I. INTRODUCTION 

Effective mass approximations have been central to an¬ 
alyze and understand band structures of materials near 
critical points in the Brillouin zone (BZ) and their ma¬ 
jor physical consequences. However, some basic formu¬ 
lae of that formalism have been misused for energy band 
dispersions that are not twice-differentiable. This has 
been further confused with band non-parabolicity effects 
in Taylor expansions. In a previous paper,— a more rig¬ 
orous theory for dealing with a broad class of energy 
band dispersions has been developed. Here we advance 
that theory by applying its formalism to properly deter¬ 
mine the density of states (DOS) and the DOS effective 
mass. The development of this formalism is necessary for 
energy dispersions at critical points that are “warped,” 
hence, conventional formulae involving second-order dif¬ 
ferentials are invalid. 

We begin by reviewing in Sec. |H] basic results of our 
previous treatment. In Sec. Ell we derive general expres¬ 
sions for the DOS and the DOS effective mass. We re¬ 
cover standard results for twice-differentiable ellipsoidal 
and hyperbolic energy dispersions in Sec. IIVI Then we 
move on to study warped energy bands that are not 
twice-differentiable. In Sec. we apply our results to 
typical degenerate valence bands in cubic materials. Con¬ 
siderable differences emerge between our proper evalua¬ 
tions of the DOS effective masses and those reported in 
original papers— We discuss those differences in Sec.IVTl 
and Sec. lVHI Subsequently, we focus on two-dimensional 
physical and mathematical models where the distinction 
between effects of band warping or “corrugation” and 


band non-parabolicity can be analytically demonstrated 
(Sec. IVHII) . Some features of those examples may be 
useful in clarifying the interplay between possible band 
warping and band non-parabolicity in non-degenerate 
conduction bands of materials that exhibit remarkable 
thermoelectric properties— Finally, we draw some con¬ 
clusions and propose further inquiries in Sec. m 

II. ANGULAR DEPENDENT ENERGY 
DISPERSION 

We have previously considered^ energy dispersions 
around a point ko in a crystal BZ in the form of 

E{kr, 6, 4>) = Eo + krQi {6, (f)) 

k^CL2{0, (jf) k^CL'iiO, (f>) . ( 1 ) 

Here, kr = |k — ko| is the radial distance between a 
generic point at k in the BZ and the point of expansion 
at ko. The latter may be any point of special interest in 
the BZ, or a “critical point,” where the energy expansion 
has a null first-order differential— The polar angles 0 
and (j) refer to the polar spherical coordinates of k — ko. 

It is essential to appreciate that Eq. (HD provides a 
much more general dispersion relation than commonly 
considered functions that admit multi-dimensional Tay¬ 
lor series expansions in Cartesian coordinates. That is 
so because Eq. (HD requires only the existence of a one¬ 
dimensional Taylor series expansion in each radial direc¬ 
tion across ko. This is a much more limited requirement 
that can be reasonably expected of any physical band 
structure that allows one-dimensional transport of quasi- 
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particles in any directionji Besides ordinary quadratic 
bands, Eq. © thus includes “warped bands,” which are 
not twice-differentiable at isolated points, by definition. 
Typical examples of warped bands derive from original 
models of Dresselhaus et al. and Kanei^d^ 

Mathematically, band warping must be unambiguously 
distinguished from band non-parabolicity. The latter de¬ 
rives from higher-order terms amid,4’) with m > 2 in 
Eq. ([T]). Conversely, band warping depends exclusively 
on the shape of the 02 ( 0 , term, which provides a di¬ 
mensionless angular effective mass surface in Rydberg 
atomic units 

For an ordinary quadratic band, associated with a 
second-order differential and its curvature, 02 ( 6 *, </>) as¬ 
sumes the form 

sin^ (9) (rrixsiii^ {(I)) + my cos'^ {(/))) cos^(9) 

^2(^5 0 ) ITJ ITY 

m^mylTle TTlzUle 

,, (2) 

in a coordinate system of principal axes, with diagonal 
effective masses , my, , while me is the ordinary 

electron mass. Any other form of 02 ( 0 , iji) corresponds to 
a warped band, which cannot be exclusively described in 
terms of diagonal effective masses. 

One may formally derive expressions for the DOS cor¬ 
responding to the general energy expansion in Eq. ©• 
In this paper, we focus on explicit DOS expressions 
for band warping, although we generalize our consid¬ 
erations at least to one type of band non-parabolicity, 
namely that of an overall energy dispersion of the form 
E{kr,9,(j)) = R{kr) f {6, (f)), where R{kr) is a monotoni- 
cally increasing function and all ai{0,(j}) = f{9,(j)). In 
this paper, we do not further consider any linear term 
in the energy expansion, thus assuming a null first-order 
differential at a “critical point.”— 


in the fornti 


E{kr,6,(j)) = Eq 


fi^k^ 

2mfi 




(4) 


In Eq. (|4]) we imply a partition of the unit (0, (j)) sphere 
in a region TZ+ where f{9, cj)) > 0 and a region TZ- 
where f{9,4’) < 0, so that 7?.+ U ??._ = S'^. This defi¬ 
nition of f{9, (j)) must refer to a single band, which may 
or may not be degenerate with other bands at kg. Typi¬ 
cally, though not necessarily^— non-degenerate bands at 
ko are not warped, corresponding to analytic maxima, 
minima, or saddle points. Conversely, degenerate bands 
are commonly warped— 

For optical transitions, the joint density of states 
(JDOS) can be similarly considered— Both conduc¬ 
tion and valence bands can be expressed as individual 
terms having the form of Eq. Q. For the JDOS we may 
then define a joint F{d, (p) as the sum of the correspond¬ 
ing two (absorbing and emitting) /-contributions. The 
same formalism that we develop in this paper for the 
DOS thus essentially applies to the JDOS as well. 

In order to proceed with the integrations in Eq. m, 
we may first scale the Cartesian coordinates by letting 
fc' = The energy dispersion thus becomes 


/2mr 


E(kr,9,cP)=Eo + k'/f(9,p). 


(5) 


We may then introduce a new variable E' = E(k') — Eq, 
so that 

^ I S{E'-{E-Eo))d^k'. ( 6 ) 
In polar coordinates we have 


III. CALCULATIONS OF DOS 


k' 



sin 9 cos 4>, 


In a crystal, the single-band DOS at energy E, within 
dE, is basically defined as 

9{E) = gsj^ J SiEik) - E) d?k, (3) 

where gs is a possible spin degeneracy, V is the volume 
of the direct-lattice primitive cell, and i5(k) represents a 
single energy band in the BZ over which the d^k integra¬ 
tion is performed. The integral may be evaluated after 
performing a transformation to (E, 9, p) coordinates. 
The delta function can thus be handled by reducing the 
integral to the surface having given energy E inside the 


k' — 

Ky - 


E' 




sm a sm < 


k' = 


E' 




cos 9. 


(7) 


Accordingly, we may perform a change of variables to 
spherical coordinates (fc(., 0 , p), where k[. is defined im¬ 
plicitly through Eq. ®. Notice that regions of positive 
E' [E > Eq) correspond to and regions of negative 
E' {E < Eq) correspond to TZ-, so that all variables in 
Eq. (I?!) are real. 

The Jacobian for this transformation is 


A. The DOS of Warped and Non-Warped Bands 

Approaching a critical point kp in the BZ, let us ig¬ 
nore band non-parabolicity effects for the moment and 
consider an energy dispersion (without any linear term) 


JiE',9,P) 


E' 


sin 0 


f{9,P) 2 /( 0 ,# 


> 0 , 


^{9,cP)&S\ 


( 8 ) 

With the change of variables to (A', 0, <()), the DOS inte¬ 
gral thus becomes 
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9{E) = g. 



E' 


sinfl 


/( 0 ,<^) 2 /( 0 , <^)2 

I- 


5{E' - {E - Eq)) d(j)d0dE'. 


(9) 


Now, if surface integrals over the unit sphere converge, 
we may split those integrals over regions of positive and 
negative /(0, (^), so that the energy integration immedi¬ 
ately yields 


E>Eo, 

\9sj^{^f"VW^[C-], E<Eo, 


where we have defined 
C'± = 


sin0 


n± 


2(±/(0,0))3/2 


dOdcj). 


( 10 ) 

( 11 ) 


However, integration over the angular variables may 
not formally converge, as in the ideal case of a saddle 
point dispersion extending to infinity*^ That is a theo¬ 
retical extrapolation, however, because the BZ is actually 
finite, and so must be any band structure within it. In¬ 
troduction of an energy-dependent cutoff parameter may 
thus be required, which should further take into account 
the onset of any significant band non-parabolicity. In 
any such case, the energy integration must be taken last, 
since C± also become functions of energy. However, the 


presence of the delta function can still make this last 
integration over energy relatively straightforward. We 
provide an example of that in Sec. IIVBI 


B. Generalization to Monotonically Non-Parabolic 
Bands 


We can readily extend the preceding formalism to en¬ 
ergy dispersions of the form E' = R{k'^) f {6, (j)), where 
R{k'j.) is any monotonically increasing function of fc(,. The 
corresponding coordinate transformations then become 


k' = R 


-1 


E' 


k' = 

" Lm</>)J 


E' 


k' = R 


-1 


E' 




sin 9 cos 
sin 9 sin </>, 
cos 9. 


( 12 ) 


The inverse function R ^ oi R must exist and it has been 
introduced in Eq. ca). The DOS thus becomes 


g{E) = 9s 


V ( 2me 


(27r)3 


s/2 poo pTZ p2'K 


R- 


' —ooJ 0 0 


E' 


I- 


sin0 

27(^ 


;5{E' - {E - Eq)) d(j3d9dE'. 


(13) 


Equation dHl) may now be regarded as a special case of 
Eq. (fT3l) . where R~^(x) = -/x. The energy integral may 
still be relatively straightforward to perform in Eq. (USD 
on account of the delta function. 


C. Two-dimensional DOS 

For further illustration, let us consider the two- 
dimensional evaluation of the DOS, according to the ex¬ 
pression 

9{E) = gsj^ J SiEik) - E) d^k. (14) 

Close to a critical point kg in the BZ, and ignoring band 
non-parabolicity, the two-dimensional energy dispersion 


becomes 


E{kr,9) = Eo + ^f{9), 

2mfi 


(15) 


where f{9) is now a function of a single angular variable. 
The Jacobian of the transformation is simply 


J 2 = 


1 


2|/(^)l’ 


(16) 


and the DOS thus becomes 
V /2me 

= “'■(sy (“fiv' 

^00 p 27 r 


Again, we must integrate over regions of positive and 
negative f(9) separately. Namely, the interval (0,27r) 
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must be split into and TZ- regions, where f{9) is 
either positive or negative, respectively. Assuming that 
corresponding 0-integrals converge, this yields 


9{E) = 


E>Eo 
E<Eo, 


where 


(18) 


C± = 


n± 


2|/(0)l 


dO. 


(19) 


IV. TWICE-DIFFERENTIABLE 
THREE-DIMENSIONAL ENERGY DISPERSIONS 

A. Ellipsoidal Energy Dispersions 


Let us begin by considering the basic case of an ellip¬ 
soidal energy dispersion of the form 








E{kx, ky, kz) — — - k^ -|- -- ky + - - k^, + Eq, ( 20 ) 


2m 


2my ^ 2mz 


where all three principal masses mx, my, and have 
the same sign. In a four-dimensional space where 
E(kx, ky, kz) refers to the fourth dimension, that provides 
a paraboloid with a minimum (maximum) if all three 
masses are positive (negative). For positive masses, we 
may rescale the corresponding coordinates as fc' = . 

The DOS is then expressed as 

V 2^/^ (mxmymz)^^^ 
g(E) = ^^ -X 




■ (27r)3 

J 5{E' -{E- Eo)) d^k' 


( 21 ) 


where E' = k^, and /(0, (p) = 1. Integration is straight¬ 
forward, yielding 


= 2tt, C- = 0. 


We thus recover the well-known result 


( 22 ) 


n(E) = < ^®( 27 r )2 R3 Eq, E > Eq, 

1 0, E< Eq. 


( 23 ) 


In the special case where all masses equal a single effective 
mass, mx = my = mz = to*, corresponding to a spherical 
energy dispersion, we have 



n ^ 
0 , 


( 2m, \3/2 


y/E — Eq, 


E > Eq 
E < Eq. 


(24) 


Comparing Equations (|23)) and (l24l) . the DOS effective 
mass is ususally defined as to* = (mxmymz)^^^. 

If all three masses are initially negative, essentially 
equivalent results can be obtained for the DOS by switch¬ 
ing the signs of all three masses to positive, while corre¬ 
spondingly switching the signs of all energies to negative 
in the preceding equations, starting with Eq. O- 
B. Hyperbolic Energy Dispersions 


A hyperbolic energy dispersion still has the form of 
Eq. (l20)l . but some of the masses have opposite signs. For 
the sake of simplicity, although without any major loss 
of generality, let us posit that mx = my = —mz = me in 
Eq. (l20l) . In that case, the angular effective mass surface 
f{9, <j)) in Eq. ([2]) and Eq. (|1]), and its corresponding 'R.± 
regions, are 


f{0,4>) = 
7^+ = 

7^_ = 


— cos(20). 


{e,cp) 

{OA) 


TT 3tT 

T < ^ < TT 
4 4 


O<0<--or^<0<7r 
4 4 


(25) 

(26) 


As expected for an endless hyperbolic dispersion, the 
surface integrals do not converge and a ‘spherical’ cut¬ 
off radius Rc must be introducedJ^ In turn, this intro¬ 
duces an energy-angle relation at Re, namely EE' = 
—ReCOs2{j ± e), where e represents the angular incre¬ 
ment from j corresponding to the intersection of the 
sphere of radius Rc and the hyperboloid of constant E'. 
Now the surface integrals based on Eq. m can be for¬ 
mally performed, yielding 
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C+{E']Rc) = 2-27r 


sin 9 


r+<= 


2(— cos 20)3/2 


d9 


— COS0 


2\/— cos 20 

CT-' 

C-{E'-R^) = 2-2tt 

COS0 




27r i?c TT y/W 
+ ■ 


{E'f/^ 


y/2y/W V2 Rc 4^2 


sin0 


0 2(cos 20)3/2 


2-\/cos20 


00 


27r i?, _ ^ ^ (-i;')3/2 

—^ . — ZTT -I- 

I- 


i?c 4^2 i?2 


(27a) 


(27b) 


We may thus perform in Eq. (|9]) the energy-dependent plete the energy integration via the delta function, thus 
angular surface integrals over 0 and (j) of C± and com- obtaining 


J 


9{E) 


9s 

V 

('2nie 

.3/2 

271- n 

h 

1 

o 

77 {E 

-Eof 1 


(27r)^ 

V 


+2^= 

' +2 Rc 

4+2 


5 

9s 

y 

1 2nie, 

.3/2 

2ir n 



Eq — E tt 

(Eo-Ef , 

(27r)^ 

(iw 

/ 

+2^= 


+2 

Rc 4+2 R'i ' " ' 


E > Eq, 
E < Eq. 


(28) 


This agrees with the form of the DOS around an Mi 
saddle point, as derived on p. 157 of Ref, fl^ for example. 

Beyond any such treatment, which is limited to twice- 
differentiable energy dispersions, we must now proceed 
to apply our angular effective mass formalism to more 
general calculations of the DOS in complex situations of 
warped energy bands. 


V. THE DOS EFFECTIVE MASS FOR WARPED 
ENERGY BANDS 


A. The DOS Effective Mass 


Since the form of Eq. (|3]) is designed to capture a band- 
warped energy dispersion at a critical point in the BZ, we 
may still use the expression in Eq. (I24p to define the DOS 
effective mass for a warped energy band minimum, or its 
straightforward modification for an energy band maxi¬ 
mum. Comparing Eq. (ITU)) with Eq. (IMl) . and recalling 
the definitions of the numerical factors given in Eq. (HU, 
we may generally define the DOS effective mass as 

Z' C+ \ 

m* = ±f—j me. (29) 


B. A Band Warping Parameter 


There are multiple ways of introducing parameters 
that provide some measures of band warping. However, 
no single parameter can be expected to account entirely 
for the full angular complexity of /(0, </>). We have previ¬ 
ously introduced one measure of band warping by defin¬ 
ing a parameter— 

( {Tr[H] - (Tr[g]))2 )V2 
{Tr[H]) 

(30) 


For the sake of illustration, let us return to a two- 
dimensional energy dispersion as in Eq. m, namely, 
E = and let us further assume that /(0) is 

positive everywhere. Then we have^ 


(•) 




00 , 


rr[iJ] = 2/(0) +2/(0+|). 


(31) 


This definition of the band warping parameter, w, essen¬ 
tially measures the coefficient of variation of the sum of 
the eigenvalues of the Hessian matrix El (0) formally ob¬ 
tained in each Cartesian coordinate system rotated by an 
angle 0. A twice-differentiable, i.e., a non-warped surface 
cannot have any variation of its quadratic form eigen¬ 
values. Hence, that must have w = 0. However, the 
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converse is not necessarily true. Nairrely, having w = 0 
is not sufficient to conclude that the energy dispersion 
is twice-differentiable. The same considerations also ap¬ 
ply to f{0) surfaces with positive and negative values in 
different regions of S^, as indicated in ^III Al 

C. The DOS Effective Masses for the Kittel Form 


value of m-hh, since that depends on the values of a and 
c parameters; and conversely. 

Our results for the angular version of the DOS effective 
mass are in fact consisteirt with those of an original pa¬ 
per by Lax and Mavroides^ if oire identifies their g{6, (p) 
with the precise angular effective mass surface /(0, p) in¬ 
troduced iir Ref. [l| and used in this context. 


As a basic illustration of our results, let us calculate the 
DOS effective masses for the heavy-hole (hh) and light- 
hole (Ih) bands described by what we may dub the “Kittel 
form,” originally derived in a ground-breaking paper,— as 


£^(k) = 

^ (Ae ± [B^k^ + C^iklkl + klkl + klkl )]. 

(32) 

Expressing that according to our Eq. ®, we obtain 
exactlyi 


f{e,<P)=A± (33) 

\JB'^ + sin^(0) [cos2(0) -I- cos^((p) sin^(0) sin^(()))]. 


In both expressions, the upper positive (lower negative) 
sign refers to the heavy (light) hole band dispersion. We 
use again Rydberg atomic units, where h^/2mg = 1. 

Although we may not be able to express it in a closed 
analytic form, each DOS effective mass for the Kittel 
form can be evaluated numerically using Eq. (ED and 
Eq. (1^91) . Let us further factorize the B parameter in 
front of the energy dispersion of the Kittel form or its 
angular effective mass surface. Contour plots of the cor¬ 
responding DOS heavy-hole effective mass, mhh, as func¬ 
tions of a = A/|i3| and c = C/\B\^ are shown in blue 
in Fig. [TJ Numerical values of mhh are given in units of 
th^B. Notice that f{9, p) becomes imaginary for some 
values of 9 and ^ if c exceeds a Cmax given by 


Cmax(a) 


4Va2 - 1 


(34) 


Contours of constant mhh thus appear to accumulate 
along a corresponding curve. It is not clear whether any 
mhh may be attained for values of a and c approaching 
Eq. (IMl) from below. 

We may also compute the band warping parameter, 
u>, for the heavy-hole band of the Kittel form, based on 
the analog of Eq. o to three-dimensional energy dis¬ 
persions. Contour plots of constant w are shown in red 
in Fig. [H Notice that, moving along curves of constant 
w, the DOS heavy-hole effective mass mhh increases with 
increasing a. Alternatively, moving along curves of con¬ 
stant mhh, the band warping parameter w increases with 
increasing c. Thus, perhaps surprisingly, a larger value of 
w does not necessarily imply either a larger or a smaller 


VI. EFFECTS OF BAND WARPING ON THE 
DOS AND THE DOS EFFECTIVE MASSES 

Given the somewhat unexpected results that we have 
obtained for the Kittel form, it is natural to question 
what effects or relations may generally exist between 
band warping and the DOS effective masses. For ex¬ 
ample, if we consider energy dispersions with angular 
contributions giving rise to finite C± in Eq. ED, then 
the only effect that band warping can have on the DOS 
is to modify that numerical factor in front of the square- 
root energy dependence in Eq. (nnD. Additional insight 
about the integrated contribution of f{9, p) in Eq. (fTTl) 
may be gained by using methods similar to that of a sta¬ 
tionary phase, that is, by looking for particular directions 
that may predominantly contribute to the overall DOS 
effective mass. In any case, it is clear that the DOS can 
be increased by increasing the effective mass given by 
Eq. (1^ in Eq. (ITU)) . 

For the Kittel form, our w parameter may also be used 
to indicate how far from spherical is the angular effective 
mass surface f{9,(j)). For example, in the plane (a,c) 



FIG. 1: Contour plots of constant the DOS heavy-hole effec¬ 
tive mass mhh, in units of meB, are shown as blue curves. 
For comparison, contour plots of constant absolute value of 
warping parameter w are shown as red curves. 
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of Fig. [U if we climb vertically along the positive c axis 
from some point, e.g. (—^312501,0), w increases.— Let 
us then compute the error between an approximate DOS 
effective mass, derived from the least-squares fit of the 
f{9, (j)) surface to a sphere, and the correct DOS effective 
mass, calculated from Eq. (HU. That error is plotted 
in Fig. [21 As expected, when c = 0, the relative error 

{■ -IS zero, because the angular eiiective mass 

surface is actually spherical. Flowever, as c and w in¬ 
crease, the relative error () increases up to 
almost 100%! Thus, at least for the Kittel form, we may 
say that w provides some indication of how far is the 
energy dispersion from being twice-differentiable . 



qqI^L. ----------- 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 


Absolute Value 
of warping parameter 

FIG. 2: Relative error versus band warping parameter, w, 
for the Kittel form. For specificity, we start at a point 
(-V312501 « -559,0) in Fig. [T] and then we increase c verti¬ 
cally. Evidently, the relative error of the DOS effective mass, 
derived from the least-squares fit of the angular effective mass 
surface to a spherical surface, increases monotonically with w. 


VII. RELATIONS TO THE LAX-MAVROIDES 
AND LAWAETZ DOS EFFECTIVE MASSES 

Lax and Mavroides^ originally proposed the correct 
idea of an angular effective mass, but they immediately 
contaminated it with questionable expansions meant to 
fit the Kittel form specifically. Their Eq. (8) and those 
at the beginning of their Sec. IIIA correspond to our 
Equations (HU) and dMl), in defining the DOS effective 
mass. However, not only is our treatment much more 
general than theirs, but it also applies more appropriately 
to the Kittel form, based on Eq. dMl). 

Our treatment of the DOS effective masses is also much 
more rigorous and clearer than that of Lawaetzi^ Using 
our generally correct expressions and integrating them 
numerically for the same values of parameters reported 
by Lawaetz for various materials, there are significant dif¬ 
ferences between our appropriate DOS effective masses 
and those artificially produced by Lawaetz. We show 


TABLE I: Comparison of the DOS effective masses for mate¬ 
rials reported in Table II of Ref. and those correctly derived 
from our Equations HU and (I29II . 


Crystal 

7l 

72 

73 

Lawaetz 

rrihd 

Correct 

mhd 

Lawaetz 

mid 

Correct 

mid 

c 

4.62 

-0.38 

1 . 

a 

a 

a 

a 

Si 

4.22 

0.39 

1.44 

0.53 

0.537 

0.16 

0.156 

Ge 

13.35 

4.25 

5.69 

0.35 

0.351 

0.043 

0.0423 

Sn 

-14.97 

-10.61 

-8.52 

0.29 

0.289 

-0.029 

-0.0297 

AlP 

3.47 

0.06 

1.15 

0.63 

0.615 

0.2 

0.195 

AlAs 

4.04 

0.78 

1.57 

0.76 

0.752 

0.15 

0.151 

AlSb 

4.15 

1.01 

1.75 

0.94 

0.953 

0.14 

0.141 

GaP 

4.2 

0.98 

1.66 

0.79 

0.786 

0.14 

0.143 

GaAs 

7.65 

2.41 

3.28 

0.62 

0.620 

0.074 

0.0739 

GaSb 

11.8 

4.03 

5.26 

0.49 

0.498 

0.046 

0.0468 

InP 

6.28 

2.08 

2.76 

0.85 

0.858 

0.089 

0.0887 

InAs 

19.67 

8.37 

9.29 

0.60 

0.600 

0.027 

0.0267 

InSb 

35.08 

15.64 

16.91 

0.47 

0.490 

0.015 

0.0147 

ZnS 

2.54 

0.75 

1.09 

1.76 

1.796 

0.23 

0.224 

ZnSe 

3.77 

1.24 

1.67 

1.44 

1.468 

0.149 

0.148 

ZnTe 

3.74 

1.07 

1.64 

1.27 

1.296 

0.154 

0.152 

CdTe 

5.29 

1.89 

2.46 

1.38 

1.466 

0.103 

0.102 

HgS 

-41.28 

-21 

-20.73 

2.78 

2.946 

- 0.012 

- 0.0121 

HgSe 

-25.96 

-13.69 

-13.2 

1.36 

1.341 

-0.019 

-0.0190 

HgTe 

-18.68 

-10.19 

-9.56 

1.12 

1.220 

-0.026 

-0.0261 


“Formalism invalid because 72 and 73 have opposite sign 


that in Table [H where we have used the following rela¬ 
tions between the A, B, and C parameters of the Kittel 
form and the 71 , 72 , and 73 parameters introduced by 
Luttinger^ 

A(7i) = - 71, 

B{i 2 ) = 272 , (35) 

C(72,73) = \/i2(7| - iD- 

Figure [3] shows the error of the DOS heavy-hole effec¬ 
tive mass estimated by Lawaetz and its correlation with 



-0.16 

t/5 

(U 

-0.18 

'c 

-0.20 


-0.22 

<L) 

B 

-0.24 

3 

-0.26 


-0.28 


0.00 0.01 0.02 0.03 0.04 0.05 
Warping parameter 
(dimensionless) 


FIG. 3: Relative error of the DOS heavy-hole effective masses 
mhd estimated by Lawaeta^ and reported in column 5 of Table 
I, as compared to our correct values, computed from Eq. HU 
and Eq. ([^ and reported in column 6 of Table I, versus the 
band warping parameter w. 























our warping parameter w for that band in various mate¬ 
rials. That error is partly the result of inconsistent se¬ 
ries expansions and truncations in procedures elaborated 
by Lax, Mavroides and Lawaetz4“— Roughly, the larger 
is warping or w, the greater is the discrepancy between 
Lawaetz’s estimate and our precise determination of the 
DOS effective mass. That error can be quantitatively 
as large as 28%. More importantly, the original lack of 
a precise definition and treatment of warped bands has 
been responsible for a lack of consistency among many 
subsequent papers and ad hoc estimates of the DOS ef¬ 
fective masses. 

To illustrate more subtle effects of band warping on 
the DOS, we shall further consider some two-dimensional 
cases where we can quantitatively control parameters 
that provide different measures of band warping, namely, 
either the w parameter that we have already introduced, 
or an alternative band warping parameter to which we 
may refer more generally as band “corrugation.” 

VIII. TWO-DIMENSIONAL CASES 


As a first case, consider a two-dimensional version of 
the Kittel form determined by setting = 0 in Eq. (15^ , 
namely, 


E{k,,ky) = ± 

= \E\^Jak^ ± ^k^ + c-^klkl). (36) 


Equivalently, by setting 0 = 7r/2 in Eq. (1551) , and then re¬ 
labeling the azimuthal angle (j) with the two-dimensional 
polar angle 0, we obtain exactly 


A. Two-dimensional Kittel form 



c=0 


c=0.4 


c=0.6 


c-0.8 


E{k,0) = 


2me 


f{0) 


\B\— -(a ± \/l -I- c2 cos2 6 sin^ 0). (37) 

2mp 


Angular effective mass planar contours of f{0) are shown 
in Eig.|5]for a given value of a and four increasing values 
of the c parameter. 


FIG. 4: Angular effective mass contours of f{9) for the two- 
dimensional Kittel form that has kz = 0, for parameter values 
of a = —1.1 and c = 0.0, 0.4, 0.6, and 0.8. The k dependence 
is exactly parabolic in every radial direction. 


In this two-dimensional case, the band warping param¬ 
eter, w, and the DOS effective mass, m*, can be expressed 
analytically, for any c < Cmax, as 


V2 


y ^2 (^2 + 8) - SA (-f )" - SV^TlE (^) E (-f) - 2 (c2 + 4) A ^ 


4™ + 4E; (-f) + 2 VH^e; 

- 1) (4a2 - c2 - 4) (^(l - a^) K (-4) + 


4(a^-l) 


+ na (a^ — l) 


ttIRI 


( a 2 _ 1 ) 3/2 ^^^2 _ c 2 - 4 . 


(38a) 


TO* = — 


(38b) 
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In Eq. (1551) . E{m), K{m), and n(n,m) denote the com¬ 
plete elliptic integral, the complete elliptic integral of the 
first kind, and the complete elliptic integral of the third 
kind, respectively, and m = sin^ a and n are their stan¬ 
dard arguments. 

Contours of constant DOS heavy-hole effective mass 
TOhh and contours of constant absolute value of warp¬ 
ing parameter w for this two-dimensional Kittel form are 
qualitatively similar to those of the full three-dimensional 
Kittel form, which was shown in Fig. [TJ In two dimen¬ 
sions w attains a maximum magnitude whenever a and 
c approach the limit of Cmax- In two dimensions, that is 

w^max = ^- « -0.8447. (39) 

We did not investigate a corresponding effect in the full 
Kittel form but we expect similar results. 


which 1/(0) I becomes increasingly smaller with increas¬ 
ing b values approaching 1~ . In fact, analytic derivations 
yield 



m* 


1 

ay/1 — 6^ ’ 


(41) 

(42) 


where m* = C+j-K. So, the band warping parameter, w, 
and the DOS effective mass, m,, are independent of each 
other, since only m* depends on a, as expected. This 
behavior may have not been anticipated, but of course 
we constructed this illustration for that purpose. 


C. Example of m? (cos^” nO -|- sin^" n@) 


B. Example of a(l-|-6cos40 ) 


ky 



b=0.1 

b=0.2 

b=0.5 

b=0.9 



FIG. 6: Angular effective mass contours of f{6) for 
a two-dimensional dispersion relation of the form E — 


^^n^(cos^"' n9 + sin^" n6) for n — 1, 2, 3, and 4. 


FIG. 5: Angular effective mass contours of f{9) for a two- 
dimensional dispersion relation of the form E — + 

bcos49), where we set a = 1, and b = 0.1, 0.2, 0.5, and 0.9. 

As a second example, we consider the two-dimensional 
energy dispersion 

fc2jL2 t2n2 

Unless 6 = 0, this function is not twice-differentiable at 
the origin exclusively, as an isolated point. In Fig. [5l 
its angular effective mass /(0) is plotted for a = 1 and 
four increasing values of 6. Since the integral in Eq. dUD 
involves |/ ^(0)|, we can indefinitely decrease m* in 
Eq. (l29l) by letting a become as small as we need. On the 
other hand, for any given value of a, we expect substan¬ 
tial contributions to m* from diagonal directions, along 


Let us now provide a more complex example where w 
steadily increases with what we may dub band “corruga¬ 
tion,” whereas rrit^ at first decreases, but then increases 
with that “corrugation.” Consider an energy dispersion 
of the form 

fc2n2 t2jL2 

E = -—fie) = -(n^ (cos^'^ nO + sin^" nO)). (43) 

2me 2me 

Again, unless n = 1, this function is not twice- 

differentiable at the origin exclusively, as an isolated 
point. In Fig. [5] we show plots of its angular effective 
mass f{9) for n = 1, 2, 3, and 4. 

The impression conveyed by Fig. [S] is that the energy 
dispersion ought to deviate more and more from being 
twice-differentiable with increasing n. We are thus led 
to regard n as a separate parameter, independent of tc, 
that may provide an alternative, albeit qualitative, mea¬ 
sure of band warping. So, we associate with n a name 
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and a notion of band “corrugation,” although that can 
hardly provide or lead to any more rigorous or general 
definition. In any case, the basic idea of “corrugation” is 
that it increases with increasing number of radial “val¬ 
leys.” One may thus expect that the DOS effective mass 
also increases correspondingly. Although that is often the 
case, it may not always be so, as we demonstrate with 
this example. In fact, we could provide many more exam¬ 
ples where the DOS effective mass increase or decrease 
with, or remains independent of, “corrugation.” 

In this example, we can still derive analytic expres¬ 
sions for the band warping parameter, w, and the DOS 
effective mass, m*. However, those expressions are fairly 
elaborate and we omit them for the sake of conciseness. 
Suffice it to say that w increases monotonically with n, 
whereas m* at first decreases with n, but then it reaches 
a minimum, after which m* increases monotonically with 
n. Corresponding plots of w and to* are shown in Fig. [T] 
This example thus demonstrates that w and to* do not 
necessarily correlate with each other, nor with the notion 
of band “corrugation.” 
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FIG. 7: Warping parameter w (in blue) and the DOS effective 
mass m* (in red and in units of me) corresponding to Eq. (IMJ. 
While w increases monotonically with n, m* decreases at first, 
but subsequently increases with n. 


or n. Consider the energy dispersion E = ^^/(0) = 

^^((n^ — lOn -I- 30) (2 -|- cos2(2n — 1)0)), whose angu¬ 
lar effective mass f{6) is plotted in Fig. [8]for n = 1, 2, 
3, and 4. We can prove that the band warping parame¬ 
ter is w = 0, independently of corrugation or n, but the 
DOS effective mass to* at first increases with corrugation, 
then it reaches a maximum at n = 3, and subsequently 
decreases monotonically for all n > 3. 

ky 

lOr 
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FIG. 8: Angular effective mass contours of f{6) for a two- 
dimensional dispersion relation of the form E = 2 m ~ 

2 m A 30) (2 + cos 2(2n — 1)0)) for n = 1, 2, 3, and 


This example demonstrates in particular that a func¬ 
tion that is not twice-differentiable can still have tc = 0. 
This prompts us to introduce in an Appendix a more re¬ 
fined definition of a band warping parameter W that cap¬ 
tures at least that type of non-differentiability. The W 
alternative to w is more elaborate but possibly more help¬ 
ful in identifying energy dispersions that are not twice- 
differentiable. However, since second-order differentiabil¬ 
ity is inherently based on multi-dimensional limits, there 
can be no single parameter whose vanishing is sufficient 
to guarantee that any particular energy dispersion is cer¬ 
tainly twice-differentiable at a point. 


D. Corrugated example with w = 0 

We have already demonstrated that the warping pa¬ 
rameter w may not necessarily increase or correlate with 
an increasing DOS effective mass nit,. In fact, looking 
back at Fig. [l] we can easily draw parametrized curves 
where w decreases while to* increases. We can also draw 
curves in Fig. [T] where w stays constant while to* either 
increases or decreases. 

Let us then provide a conclusive two-dimensional ex¬ 
ample that has 1/3 = 0, although the energy disper¬ 
sion is not twice-differentiable, and rrif still decreases 
at first, and then increases with increasing corrugation 


IX. CONCLUSIONS 

We have applied the angular effective mass formalism 
introduced in Ref. [l| to study the density of states in 
warped and non-warped energy bands at critical points 
in the Brillouin zone. First we have verified ordinary 
results for ellipsoidal and hyperbolic energy dispersions. 
Then we have generalized the expression of the DOS to 
account for general band warping and monotonically in¬ 
creasing non-parabolic energy dispersions. Band warp¬ 
ing may or may not increase the DOS effective mass. An 
intuitive notion of greater band “corrugation,” referring 
to energy dispersions that deviate “more severely” from 
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being twice-differentiable at an isolated critical point, 
may also vary independently of the corresponding DOS 
effective mass and band warping parameter. We have 
demonstrated these effects through investigation of va¬ 
lence band energy dispersions in cubic materials, showing 
the superiority of the angular effective mass treatment of 
the DOS effective masses compared to that of original 
papersi^^— 

We have further considered certain two-dimensional 
physical and mathematical examples that may be rele¬ 
vant to studies of band warping in heterostructuresi^^— 
and surfaces^!. These examples may also be useful in 
clarifying the interplay between possible band warping 


and band non-parabolicity for non-degenerate conduc¬ 
tion band minima in thermoelectric materials of corre¬ 
sponding interest.— “— 
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X. APPENDIX 

If f{x, y) is a twice-differentiable function of two Carte¬ 
sian variables at a point a, then the second-order direc¬ 
tional derivative at a is defined as 


J 


Dif{ax,ay) = Imr 


f{ax + 2tcos{9),ay -I- 2tsin(0)) - 2f{ax + tcos{9),ay -I- tsin(0)) J- f{ax,ay) 

■ 


(44) 


This directional derivative can also be expressed as a 
linear combinations of second-order partial derivatives 
along the coordinate x- and y- axes at a. 

If A is an orthogonal matrix with determinant -1-1, 
whose first column is derived from the x-axis rotated into 
a new direction by an angle 9, and x is the unit vector 


along the original x-axis, then one can show that 

Djfiax, ay) = xTATiL(6»)/(a)Ax, (45) 

where H{9) is the Hessian matrix of ordinary second- 
order partial derivatives at a. 
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If /(x, y) is twice-differentiable at a, the band warp¬ 
ing parameter that we have previously introduced must 
vanish.— Namely, ic = 0 is a necessary condition for 
second-order differentiability at a critical point. How¬ 
ever, ui = 0 is not a sufficient condition for second-order 
differentiability at a critical point. Expecting that any 
single parameter could capture the full complexity of 
f{x,y) around a would indeed be asking too much. 

For example, consider g{x,y) defined as zero at the 
3x® — 'dx^y'^ + 21x'^y'^ + y® 

origin and g{x,y) = -—-- every- 

(x2 + y2) 

where else. That has w = 0, although g{x,y) is not 
twice-differentiable at the origin, exclusively. In fact, its 
angular function is (2 -|- cos(60))r^, which corresponds to 
the last example given in the previous text for n = 1. 
More generally, any function of the form r'^f{6), where 

OO 

/(6I) = oo + ^ a„cos(2(2n - 1)0), (46) 

n—1 

shares the same peculiarity. 

These considerations prompt us to consider other mea¬ 
sures of band warping for functions that are not twice- 
differentiable. Consider, for example, the difference be¬ 


tween the second-order directional derivative at a and the 
correspondingly rotated Hessian matrix element, namely, 

S{0) = Dgf{ax, Oy) - (47) 

According to Eq. ([45]), if f{x,y) is twice-differentiable 
at a, then 6(9) must vanish for any 9. Unfortunately, 
the converse still cannot provide a sufficient condition in 
general. Nevertheless, detecting a non-vanishing S{9) at 
any 9 may provide a more refined tool to discover whether 
f{x,y) is not twice-differentiable. That is accomplished 
by evaluating the single parameter 

/ 1 \ 

^=[■^1 ■ (48) 

For the previous example of g(x, y) with w = 0, Eq. ((4^ 
indeed provides a non-vanishing W = 34. Thus, if we 
use W rather than w, we can conclude that g(x, y) is not 
twice-differentiable at the origin. Most generally, how¬ 
ever, even W = 0 cannot guarantee that any particu¬ 
lar function f(x,y) is certainly twice-differentiable at a 
point. 



